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基于 自 适 应 双 国 值 的 计 步 算法 


王 网 ， 茧 ， 敏 "， 周 清 峰 
(合肥 工业 大 学 计算 机 与 信息 学 院 , 合肥 230009) 


摘 要 : 针对 现 有 计 步 算法 对 不 同 运动 状态 适应 性 差 的 问题 ， 设 计 了 一 种 基于 自 适 应 阔 值 的 计 步 算法 。 该 算法 首先 
通过 智能 手 环 的 内 置 三 轴 加 速度 传感器 采集 用 户 在 慢 走 、 快 走 以 及 跑步 三 种 步行 频率 的 加 速度 数据 ， 经 过 五 点 滤波 
预 处 理 ， 在 自 适应 时 间 窗 内 检测 波峰 波 谷 ， 再 将 波峰 均值 和 波 谷 均值 的 平均 值 作 为 上 闪 值 ， 波 谷 均 值 作为 下 冰 值 ， 
通过 动态 阅 值 判定 步 数 ， 最 后 根据 行走 振幅 和 行走 频率 存在 的 规律 性 进行 假 步 检 测 。 测 试 表明 ， 该 算法 对 于 不 同 的 
用 户 在 三 种 不 同 频率 下 平均 计 步 精度 可 以 达到 91.88% 以 上 。 
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Pedometer method based on self-adaptive double threshold algorithm 


Wang Lan, Peng Mini, Zhou Qingfeng 
(School of Computer Science & Information Engineering, Hefei University of Technology, Hefei 230009, China) 


Abstract: To solve the problem that the existing step counting algorithms have poor adaptability to different motion states, 
this paper proposed a step counting algorithm based on self-adaptive threshold. By using the 3-axis accelerometer sensor in 
smart bracelet, the proposed algorithm first performed the acceleration data collection when users walked at three 
frequencies including slow walking, fast walking and running. After five-point filter preprocessing, this algorithm detected 
peaks and valleys in the self-adaptive time window. Then it taken the average of the peak mean and valley mean as the 
upper threshold, and used the valley mean as the lower threshold. Based on this, the algorithm adopted a dynamic threshold 
analysis to achieve the step counting. Finally, in this algorithm, false steps could also be detected according to the regularity 
of walking amplitude and frequency. The test shows that the achievable average step counting accuracy of the proposed 
algorithm is above 91.88% for different users at three different walking frequencies. 

Key words: self-adaptive; accelerometer sensor; dynamic threshold; step counting accuracy 


je 于 每 个 人 的 运动 状态 不 同 ， 步 幅 的 差别 会 比较 大 ， 所 以 该 算 
法 存在 的 缺陷 在 于 对 不 同 运动 状态 没有 普 适 性 。 虽 然 基于 阔 
近年 来 ， 随 着 人 们 越 来 越 注 重 身体 健康 问题 ， 健 身 日 益 。“ 值 的 计 步 方法 适用 性 是 很 强 的， 但 是 对 目前 而 言 ， 还 是 很 难 
受到 关注 。 由 于 步 检测 和 计数 可 以 为 个 体 的 身体 活动 水 平 提 ”确定 一 个 统一 的 阐 值 来 适应 不 同 的 用 户 和 多 样 的 行走 状态 
供 基本 的 衡量 I， 使 得 利用 智能 手 环 等 可 穿戴 式 智能 设备 来 。 nol。 
“三 。 进行 计 步 和 获取 其 他 生理 数据 成 为 了 研究 热点 。 目 前 ， 智 能 十 对 现 有 算法 对 不 同 运 动 状态 适应 性 差 和 无 法 适应 不 同 
8 手 环 的 种 类 繁多 ， 计 步 方法 也 多 种 多 样 。 当 前 计 步 算法 均 是 “用户 行走 特征 的 问题 ， 设 计 了 一 种 基于 自 适 应 双 阔 值 的 计 步 
基于 加 速度 传感器 展开 的 。 文 献 D] 提 出 了 自 相关 分 析 法 ， 通 “算法 。 该 算法 先 将 采集 到 的 加 速度 数据 进行 符号 统一 化 ， 经 
过 计算 自 相关 系数 和 标准 差 ， 设 置 经 验 阔 值 来 判断 空闲 和 行 。 过 五 点 滤波 法 对 加 速度 数据 进行 预 处 理 操作 ， 接 着 在 自 适应 
走 的 状态 ， 但 是 该 算法 未 统计 当 人 变速 运动 导致 每 步 周期 发 ”时间 窗 内 根据 波峰 波 谷 的 值 分 别 计算 上 效 值 和 下 阔 值 ， 再 通 
生变 化 时 的 计 步 效果 。 文献 3] 提出 一 种 基于 加 分 减 分 策略 的 “， 过 动态 阔 值 来 判定 步 数 ， 最 后 根据 行走 振幅 和 行走 频率 存在 
计 步 算法 ,在 固定 时 间 窗 内 比较 相 邻 两 个 采样 点 幅 值 的 大 小 ， ”的 规律 性 进行 假 步 检 测 , 可 以 准确 地 排除 非 行 走 因素 (起 立 、 
通过 比较 得 到 的 关系 来 统计 分 数 ， 当 分 数 达到 一 定 要 求 ， 则 。” 坐 下 、 打 字 、 训 桌子 等 ) 对 正常 步 数 的 干扰 ， 间 在 提高 计 步 算 
计 为 1 步 ， 但 当 人 物 步 频 变化 时 , 采用 固定 时 间 窗 往往 不 能 。 ”法 对 不 同 运动 状态 以 及 不 同 用 户 行走 特征 的 适应 性 。 实 验 结 
精确 计 步 四 。 文 献 [5] 介 绍 了 一 种 自 适应 的 波峰 检测 算法 ， 在 。” 果 表 明 ， 本 算法 在 不 同 用 户 的 惕 走 、 快 走 、 跑 步 等 状态 下 都 
正常 态 和 非 正 常态 时 计 步 准确 度 可 以 达到 99% 以 上 ， 但 是 并 ”能 获得 较 高 的 计 步 精度 。 
未 统计 算法 对 不 同 用 户 行走 特征 的 适应 性 。 根 据 阔 值 判定 来 3 RE 
进行 计 步 是 目前 常见 的 一 种 计 步 方法 ， 通 过 给 用 户 系 上 一 个 1 ” 人 体 行走 模型 与 步 数 检测 算 ; 
惯性 传感器 ， 当 测 得 的 加 速度 传感器 数据 超过 闵 值 ， 用 户 的 人 体 在 行走 过 程 中 ， 手 辟 在 一 定 方向 上 会 有 规律 性 的 摆 
步 数 将 会 增加 一 步 ka。 文 献 [g] 提 出 先 求 出 动态 峰值 和 谷 值 。 动 ， 本 文 可 以 把 这 种 摆动 看 做 是 一 种 钟 摆 运 动 1。 如 图 1 所 
> 差 ， 再 将 结果 和 固定 阔 值 比较 来 判断 步 数 是 否 增加 。 但 示 , 手臂 从 A 到 卫 就 表示 一 个 完整 的 步伐 ， 如 此 双手 会 有 
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个 不 断交 替 的 运动 。 
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1.1 行走 参数 分 析 


图 1 “手臂 摆动 模型 


Fig.1 Arm swing model 


2 ”基于 自 适 应 双 阅 值 的 计 步 算法 设计 


上 述 算法 的 不 足 之 处 在 于 人 体 在 行走 过 程 中 不 可 能 一 直 
以 一 个 速度 前 进 ， 因 此 每 一 步 的 周期 不 尽 相 同 。 当 人 行走 速 
度 发 生变 化 时 ， 算 法 的 性 能 明显 下 降 ， 因 此 提出 了 基于 自 适 
应 双 阔 值 的 计 步 算法 。 考 虑 算法 对 不 同 用 户 行走 特征 的 适应 
性 ， 采 集 了 多 个 用 户 的 行走 数据 。 算 法 的 基本 流程 是 数据 采 
集 、 数 据 预 处 理 、 计 算 自 适应 时 间 窗 长 度 ， 利 用 计算 所 得 的 
时 间 窗 长 度 进行 步 数 检测 以 及 假 步 检测 。 
2.1 数据 采集 
本 算法 的 数据 是 由 内 置 在 智能 手 环 内 的 三 轴 加 速度 传 感 
器 采集 得 到 的 。 这些 加 速度 数据 可 以 在 安 卓 手 机 APP 的 控制 


从 行走 参数 角度 ， 行 人 肢体 上 
数 如 加 速度 、 角 速度 和 旋转 角速度 
人 在 正常 行走 时 进入 


期 性 地 运动 会 导致 运动 参 
等 发 生 周期 性 的 变化 。 行 
， 在 此 过 程 中 ， 切 向 加 速 


ES 


手臂 摆动 阶 


度 由 小 变 大 ， 再 由 大 变 小 。 因 此 ， 通 过 分 析 人 体 行走 过 程 中 


1.2 步 数 检 测算 法 


波峰 检测 算法 。 
1.2.1 自 相 关 分 析 法 
自 相 关 分 析 法 


常见 的 步 数 检 测 


的 加 速度 数据 的 变化 可 以 进行 步 数 检测 [23。 


算法 有 自 相 关 分 析 法 、 加 减 分 计 步 算法 、 


是 利用 人 体 在 连续 运动 过 程 中 产生 的 加 速 


度数 据 序列 之 间 的 相似 性 来 进行 步 数 的 计算 所 。 人 体 在 行走 
过 程 中 ， 加 速度 数据 存在 一 定 的 周期 性 。 自 相关 分 析 法 就 是 
通过 计算 加 速度 数据 在 前 后 两 个 周期 中 的 相关 性 来 计 步 的 。 


计算 整体 加 速度 后 利用 式 (1) 求 整体 加 速度 的 自 相关 系数 。 
> [xm+R) umt)x(x(m+tk+t) um+t,t))] 


w(m,1)= 


1x6(m,t)xo6(m+t,t) (1 


其 中 :x(mt) 和 5(mt) 分 别 表 示 的 是 当前 加 速度 序列 


{XR),Xx(K+D,…,x(k+1 一 D)} 的 均值 和 标准 差 。 当 采样 周期 t 


接近 行人 的 行走 周期 时 ， 自 相关 系数 wm) 的 值 接近 1。 然 
司 一 行人 在 不 同时 刻 的 行走 周期 是 不 一 样 的 ， 


而 不 同行 人 或 者 
所 以 t 是 一 个 变量 。 


[FHl 


i 的 动态 确定 使 用 加 框 算法 来 实现 ， 将 t 


选 定 一 个 范围 区 间 twin-tmax, 通过 式 (2) 计 算 自 相关 系数 p(mt) ， 
当 (mt) 达到 最 大 值 时 的 上 值 即 为 该 次 行走 的 周期 。 


PUmD=max2es (wl,t)) (2) 


=tmin 


1.2.2 加 减 分 计 步 算法 


该 方法 的 基本 思想 是 对 行走 过 程 中 加 速度 数据 的 峰值 进 


he 


此 算法 认为 是 ， 


行 检测 。 算 法 在 针对 
再 依靠 分 析 阅 值 的 方法 ， 而 是 通过 计 分 机 制 。 


时 间 窗 内 峰值 可 能 存在 较 多 的 情况 ， 不 


在 每 次 计 分 时 ， 首 先 比较 相 邻 加 速度 值 


大 小 。 如 果 左 侧 相 邻 


个 加 速度 数据 增 大 则 加 1 分 , 反之 减 


1 分 ; 如 果 右 侧 相 邻 两 


个 加 速度 数据 依次 减 小 则 加 1 分 , 反 


之 则 减 1 分 ;如 采 
则 认为 走 了 一 步 。 
1.2.3 波峰 检测 算法 


具体 来 说 ， 就 是 


个 时 间 窗 的 计 分 总 和 大 于 某 一 计 分 阔 值 ， 


分 析 人 体 行走 时 加 速度 数据 出 现 的 正弦 


波 的 特性 ， 可 以 通过 检测 加 速度 正弦 波 的 波峰 或 者 波 谷 的 个 
数 来 识别 步伐 ， 波 峰 检测 算法 正 是 通过 检测 加 速度 波峰 来 实 
现 计 步 的 功能 ， 若 连续 检测 到 两 个 波峰 则 记 为 一 步 0 


假设 当前 算法 


窗口 大 小 为 xn， 窗口 内 的 加 速度 序列 


402D={fxzK+LD,xK+20 x+z2D ， 则 满足 条 件 的 加 速度 峰 


值 Vn 可 以 表现 为 如 


Xpeat = XK+m,1) > xCK+i,D 
m<nandiell,m)V (m,nl] 


式 (3) 所 示 。 
GB3) 


下 通过 蓝牙 协议 传输 至 手机 中 ， 进 而 可 以 通过 WiFi 或 数据 
线 将 加 速度 数据 传送 至 电脑 上 进行 进一步 的 处 理 。 和 截取 了 阁 
分 采集 数据 ， 其 中 从 左 到 右 包 括 序号 、X 轴 加 速度 、Y 轴 加 
速度 、Z 轴 加 速度 、 采 集 时 间 ， 如 表 1 所 示 。 
表 1 部 分 采集 数据 格式 
Table 1 Partial collected data format 


~ 


序号 X 轴 Y 轴 Z 轴 年 /月 / 时 间 
58 -22 -5 16 2017/1/3 21:47:38:314 
59 -22 -6 16 2017/1/3 21:47:38:314 
60 -26 -10 17 2017/1/3 21:47:38:314 
61 -28 -10 15 2017/1/3 21:47:38:314 
62 -31 -10 14 2017/1/3 21:47:38:314 
63 -32 -9 11 2017/1/3 21:47:38:332 


2.2 数据 预 处 理 
利用 MATLAB 分 别 画 出 了 X、Y、2Z 三 个 轴 的 加 速度 数 
据 波 形 ， 如 图 2 所 示 。 为 了 方便 数据 处 理 ， 首 先 通过 求 取 合 
加 速度 的 方式 将 数据 进行 符号 统一 化 。 分 析 可 知 大 部 分 的 Y 
轴 数 据 波 形 振幅 都 较 小 ， 出 现 这 种 情况 的 原因 是 人 体 在 正常 
行走 过 程 中 手臂 一 般 在 上 下 前 后 两 个 方向 上 摆动 ， 左 右 摆动 
幅度 较 小 。 计 算 X,Z 两 个 轴 的 合 加 速度 ， 波 形 如 图 3 所 示 。 
求 取 合 加 速度 如 式 (4) 所 示 。 
Xs = VX2+Z2 (4) 
于 噪声 以 及 不 同人 之 间 行 走 频率 的 差异 ， 数 据 波形 有 
很 多 毛刺 。 为 了 消除 毛刺 对 计 步 精度 的 影响 ， 必 须要 对 数据 
进行 滤波 。 本 文 所 采用 的 滤波 方法 是 五 点 平滑 滤波 ， 依 此 得 
到 平滑 波形 ， 如 图 4 所 示 。 五 点 滤波 公式 如 下 ， 其 中 X(n) 表 
示 第 n 个 采样 点 。 滤 波 如 式 (5) 所 示 . 
ee (5) 


四 
-20 上 
E 
-40 上 
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2 三 轴 加 速度 数据 波形 


Fig.2 Data waveform of tri-axial acceleration 
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图 3 二 轴 合 加 速度 数据 波形 


Fig.3 Data waveform of biaxial acceleration 
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图 4 二 轴 合 加 速度 滤波 后 的 数据 波形 
Fig.4 Data waveform of biaxial acceleration after filtering 


2.3” 自 适应 时 间 窗 设计 
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据 的 有 效 性 ， 一 般 截 取 中 间 部 分 的 数据 。 本 算法 所 采用 的 数 
据 是 基于 第 二 种 方法 获得 的 ， 采 样 个 数 选取 为 300。 计 算 公 
式 如 式 (6) 所 示 。 

Iss > op dimg IN (6) 


2.4 步 数 检测 
2.4.1 检测 波峰 波 谷 
在 一 定时 间 窗 长 度 内 ， 对 加 速度 数据 进行 滤波 等 预 处 理 
操作 ， 以 去 除 高 频 噪 声 产生 的 毛刺 。 设 加 速度 经 过 处 理 后 为 
XGO)， 从 第 二 个 数据 起 开始 进行 波峰 波 谷 检测 。 
如 果 时 间 窗 的 加 速度 数据 大 小 满足 
X(n-N+l)-X(n-N)>0 
X(n)-X(n-1)>0 
| (7) 
X(n+N+l)-X(nt+N)>0 
即 视 X(w) 为 时 间 窗 内 的 一 个 峰值 ， 并 将 其 存储 在 max_data 数 
组 中 ; 
同 理 ， 如 果 时 间 窗 的 加 速度 数据 大 小 满足 
XO-N+D-XO-N<0 
X(n)-X(n-1)<0 
0 (8) 
Xn+N+1l)-X(n+N)<0 
即 视 X(w) 为 时 间 窗 内 的 一 个 谷 值 ， 并 将 其 存储 在 min_dara 数 
组 中 。 
2.4.2 双 国 值 
算法 采用 的 计 步 方法 是 双 阔 值 分 析 法 。 所 选取 浆 值 是 
否 合适 对 计 步 精度 有 决定 性 的 影响 。 双 阔 值 即 以 时 间 窗 内 波 
峰 均 值 与 波 谷 均值 的 平均 值 作 为 上 阔 值 ， 波 谷 均值 作为 下 阔 
值 . 如 果 采 用 波峰 均值 为 上 闭 值 就 会 有 很 大 的 几率 产生 漏 检 。 
阀 值 计算 公式 如 下 : 


Ty =(H+L)/2 (9) 
7,=L (10) 


不 同 的 人 有 不 同 的 步行 频率 ， 如 果 按 照 固 定 的 时 间 窗 长 


其 中 : H 是 时 间 窗 内 的 峰值 均值 ;L 是 时 间 窗 内 的 谷 值 均值 ; 


度 来 测量 步 数 肯定 会 产生 误差 ， 因 此 要 找 出 一 种 方法 来 适应 
不 同人 的 步行 频率 04。 
本 算法 所 采用 的 方法 就 是 自 适应 时 间 窗 .所 谓 自 适应 时 
间 窗 ， 就 是 通过 分 析 采 样 数据 计算 每 一 有 效 步 数 的 采样 个 数 
的 均值 ， 并 将 其 作为 检测 用 户 步 数 的 时 间 窗 长 度 。 在 计算 时 


7s 是 当前 时 间 窗 内 的 上 阐 值 ，Z 是 当前 时 间 窗 内 的 下 立 值 。 
不 过 由 于 时 间 窗 内 极 值 的 个 数 不 定 ， 所 以 在 计算 闵 值 的 
时 候 要 遵循 一 定 的 方法 。 
1) 多 极 值 
当时 间 窗 内 有 两 个 或 两 个 以 上 的 峰值 或 谷 值 时 ， 可 以 直 


A EE 


间 窗 长 度 的 过 程 中 ， 先 设置 一 个 时 间 窗 长 度 初 值 。 之 所 以 这 
两 个 原因 : 一 是 方便 对 采样 数据 进行 步 数 检测 ， 以 求 
有 效 步 数 的 采样 个 数 ， 二 是 由 于 样本 数据 较 少 ， 先 设 
才 间 窗 长 度 不 会 对 计 步 精度 有 较 大 影响 ， 不 过 在 计算 过 程 
能 会 存在 因为 采集 的 数据 不 当 而 导致 时 间 窗 过 大 。 为 了 
时 间 窗 长 度 过 大 而 导致 计 步 精度 降低 的 问题 ， 本 算法 在 
[多 次 数据 测试 后 定义 了 一 个 取 时 间 窗 长 度 的 条 件 ， 即 当 


+ 


接 计 算 波峰 波 谷 的 均值 , 继而 得 到 双 阔 值 , 如 式 (11)(12) 
所 示 。 


H= Smax_datali) /No (11) 
i=l 


工 = min_daraG) /Nan (12) 
i=l 


其 中 :及 是 时 间 窗 内 的 峰值 均值 ; 工 是 时 间 窗 内 的 谷 值 均值 ; 


EF 防 骨 卫 欠 寻 人 
没 汪 涂 瑟 豆 太 夸 


项 


的 自 适 应 时 间 窗 长 度 大 于 一 定 值 时 ， 就 将 时 间 窗 长 度 减 


lls < 
T 


I 
o 


算法 通过 检测 样本 数据 的 有 效 步 数 ， 利 用 双 阔 值 分 析 法 
计算 出 每 一 步 的 采样 点 数 ， 并 将 其 存放 在 数组 中 ， 最 终 经 过 
假 步 检测 后 计算 出 每 一 有 效 步 数 的 平均 采样 点 数 ， 这 个 采样 
点 数 就 是 本 文 所 计算 的 自 适 应 时 间 窗 长 度 。 设 置 的 时 间 窗 长 
度 是 经 过 多 次 数据 测试 后 根据 具体 的 计 步 精度 确定 的 ， 减 弱 
了 其 对 计 步 精度 的 影响 。 另 外 计算 时 间 窗 的 数据 可 以 通过 两 
种 方法 获得 : 一 种 是 单独 采样 ， 采 集 每 个 人 在 不 同 频率 下 行 


二 


pa 


max_datali) 是 当前 时 间 窗 内 的 波峰 值 ，min_data(?) 是 当前 时 间 
窗 内 的 波 谷 值 ， Ns 是 时 间 窗 内 波峰 的 个 数 ， wu 是 时 间 窗 
内 波 谷 的 个 数 。 
2) 单 极 值 
当时 间 窗 内 波峰 或 者 波 谷 只 有 一 个 时 ， 取 上 个 时 间 窗 的 
对 应 的 阔 值 与 本 时 间 窗 的 极 值 的 均值 作为 本 时 间 窗 内 的 阔 值 ， 
如 式 (13 (14) 所 示 。 
H=tmp_H+max_data (13) 
L=tmp_L+min_data (14) 


走 一 定 步 数 或 一 定时 间 后 的 数据 作为 样本 数据 ， 另 一 种 方法 
是 截取 一 段 进 行 计 步 的 用 户 数据 作为 样本 数据 ， 为 了 保证 数 


其 中 : tmp_H 是 上 一 个 时 间 窗 内 的 峰值 均值 ，wmp_ 工 是 上 一 
个 时 间 窗 内 的 谷 值 均 值 ，max_data 是 当前 时 间 窗 内 的 波峰 值 ; 
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min_data 是 当前 时 间 窗 内 的 波 谷 值 。 
3) 零 极 值 本 
当时 间 窗 内 没有 波峰 或 者 波 谷 ， 取 上 个 时 间 窗 的 对 应 的 . | 
阔 值 作为 本 时 间 窗 内 的 阔 值 ， 如 式 〈15)《〈16) 所 示 。 
H=tmp_H (15) 
L=tmp_L (16) } 个 
其 中 : wmp_H 是 上 一 个 时 间 窗 内 的 峰值 均值 ，wrp_L 是 上 一 a 
个 时 间 窗 内 的 谷 值 均值 。 
2.4.3 步 数 检测 30 上 伪 波 从 
步 数 检测 以 采集 到 的 经 过 处 理 的 加 速度 信号 波峰 值 大 于 
上 阔 值 Tw 且 此 峰值 下 面 的 一 到 两 个 采样 数据 小 于 下 阔 值 7 28T 
时 记 为 一 步 。 图 5 所 示 为 步 数 检测 与 假 步 检测 的 整体 流程 。 W177 
窗 内 的 峰值 1/50s 
大 于 上 阐 值 T! . 、 、 
图 6 波峰 检测 中 伪 波 峰 的 干扰 


Fig.6 Interference ofpseudo-peaks during peak detection 


32 
波形 峰值 下 面 的 值 
是 否 小 于 下 阔 值 T 31.9 
3 波峰 波 谷 差 值 小 于 1 
31.7 天 
步 数 加 1， 将 峰值 下 标 保 3161 
存在 数组 tmp_dim 中 
如 别 三 
| 
31.4 | I 
| 
SS ‖ 
31.2 
将 峰值 下 标 保存 在 数组 a 
tmp_dim 中 31 h h h h h h h h h 
0 50 10 150 200 250 30 350 400 450 500 
1/50s 
图 静止 状态 下 二 轴 合 加 速度 的 数据 波形 
步 数 减 再 波 合 的 差 人 图 7 静止 状态 下 二 轴 合 加 速度 的 数据 波形 
2 是 否 大 于 1 Fig.7 Data waveform of biaxial acceleration in static state 


3 ”实验 与 分 析 


在 生活 中 ， 每 个 人 的 步行 频率 都 不 一 样 ， 在 采集 数据 时 
同一 个 人 在 行走 过 程 中 也 不 可 能 一 直 以 相同 的 速度 前 进 。 要 
图 5 步 数 检测 与 假 步 检测 流程 使 算法 具有 实用 性 ， 就 需要 能 够 适应 各 种 步行 频率 。 为 了 达 
Fig.5 Flowchart of step number detection and false step detection 到 这 个 目的 , 在 数据 采集 过 程 中 , 用 不 同 用 户 的 慢 走 、 快 走 
2.5 假 步 检测 跑步 这 三 种 方式 来 代表 不 同 频率 。 为 了 使 算法 更 具有 代表 性 ， 
2.5.1 行走 频率 测试 了 一 些 在 特殊 情况 下 如 将 手 环 放 在 口袋 内 ( 即 不 甩 臂 慢 
于 人 的 正常 行进 (包括 慢 走 、 快走 和 跑步 ) 步 频 为 0.2~5 走 ) 的 数据 。 
Hzb5， 本 算法 的 数据 采样 频率 /是 S0 Hz， 所 以 在 步 数 检测 3.1 算法 测试 
阶段 所 检测 到 的 每 两 步 之 间 的 采样 点 数 要 在 0.2 六 5f 之 间 (/ 为 通过 对 采集 到 的 数据 分 析 发 现 ， 在 运动 状态 下 波形 幅度 
手 环 的 数据 采集 频率 )， 和 否则 视 为 假 步 ， 并 且 此 峰值 不 计 步 。 明显 比 静 止 状态 下 的 要 大 ， 因 此 在 算法 拟 合 中 拟 采用 波形 幅 
公式 如 (17) 所 示 。 度 的 大 小 来 区 分 运动 与 静止 状态 。 完 成 算法 后 ， 进 行 了 三 轴 
dim=tmp _ dim(?) -tmp_ dim(i-D) (17) ”数据 测试 。 对 X、Y、2Z 三 轴 加 速度 的 合 加 速度 进行 分 析 ， 
其 中 : 数组 tmp_dim 中 存放 着 有 效 峰 值 ， 当 两 个 峰值 之 间 的 。 ”得 到 的 测试 结果 如 表 2 所 示 。 
差 值 即 采样 点 数 dim 满足 0.2f<dim<5f， 则 将 其 暂时 计 为 有 效 对 算法 进行 改进 ， 由 于 人 在 行走 过 程 中 手臂 几乎 在 一 个 
步 数 ， 否 则 为 无 效 步 数 。 图 6 所 示 为 波峰 检测 中 伪 波 峰 干 扰 ”平面 上 摆动 (上 下 前 后 方向 ,左右 波动 较 小 )， 所 以 ,为 了 尽 
的 示意 图 。 可 能 的 排除 其 他 因素 的 干扰 , 将 数据 波形 振幅 都 较 小 的 Y 轴 
2.5.2 行走 振幅 不 加 以 考虑 ， 第 二 次 测试 只 选择 分 析 与 Z 轴 的 数据 ， 并 称 
通过 波形 发 现 ， 人 体 在 运动 与 静止 状态 下 所 采集 到 的 数 ”为 二 轴 数 据 测试 ， 由 式 (4) 求 得 合 加 速度 数据 ， 测 试 结果 如 表 
据 波形 幅度 并 不 一 样 ,其 波峰 波 谷 的 差 值 都 在 一 定 的 范围 内 。 ”3 所 示 。 
经 过 大 量 的 数据 测试 ， 将 波峰 波 谷 的 差 值 小 于 1 的 状态 定义 从 测试 结果 可 以 看 出 ， 对 于 使 用 二 轴 加 速度 实现 计 步 ， 
为 静止 状态 ， 如 图 7 所 示 。 这 种 方法 有 效 地 排除 了 用 户 在 静 ”除了 将 手 环 放 在 口袋 中 行走 时 准确 率 有 所 降低 外 ， 在 用 其 他 
止 状态 下 对 计 步 精度 的 干扰 问题 状态 行走 时 ， 准 确 率 有 明显 的 提升 ， 而 且 将 手 环 放 在 口袋 中 
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于 


岚 ， 等 : 


之 所 以 准确 率 会 有 所 降低 ， 是 因为 手 环 内 的 三 轴 加 速度 传 感 


器 方向 无 法 确定 。 在 实际 生活 中 , 一 般 手 环 都 会 带 在 手臂 上 ， 
所 以 这 种 状态 下 精度 降低 是 可 以 接受 的 。 
表 2 三 轴 数 据 测试 结果 
Table 2 Testresults of 3-axis data 
必要 自 适 应 30 
测试 项 实际 步 数 。 ，-，,，，， 误差 步 数 ( 步 ) ”准确 率 
实 显 步 数 
昌 臂 慢 走 200 202 2 99.00% 
100 88 -12 88.00% 
100 88 2 起 88.00% 
登山 慢 走 100 107 人 93.00% 
有 辟 快 走 200 154 -46 77.00% 
200 162 -38 81.00% 
不 甩 臂 慢 走 200 205 5 97.50% 
200 178 -22 89.00% 
200 179 -21 89.50% 
100 96 -4 96.00% 
跑步 100 76 -24 76.00% 
100 103 3 97.00% 
表 3 二 轴 数 据 测试 结果 
Table 3 Test results of 2-axis data 
测试 项 目 。 实际 步 数 。 、 详 39 误差 步 数 ( 步 》 。 准确 率 
实 显 步 数 
甩 臂 慢 走 200 182 -18 91.00% 
100 96 -4 96.00% 
100 92 -8 92.00% 
登山 慢 走 100 99 | 99.00% 
电 辟 快走 200 186 -14 93.00% 
200 180 -20 90.00% 
不 甩 辟 慢 走 200 186 -14 93.00% 
200 185 -15 92.50% 
200 173 -27 86.50% 
100 89 -11 89.00% 
跑步 100 80 -20 80.00% 
100 94 -6 94.00% 


3.2 对比 分析 


算法 完成 后 ， 将 算法 与 小 米 手 环 2 


的 精度 作 了 对 比 。 对 


比方 法 是 将 小 米 手 环 2 和 A10 手 环 (算法 经 过 10 次 修改 与 调 


试 ， 故 称 为 A10 手 环 ) 分 别 戴 在 左右 手 上 ， 在 走 过 
数 后 记录 下 小 米 手 环 2 的 步 数 ， 
的 数据 经 过 算法 计 步 得 到 的 步 数 ， 然 后 将 P 


比 数据 如 表 4 所 示 。 
从 对 比 结果 可 以 看 
度 比 小 米 手 环 2 的 精度 高 ,不 


一 定 的 步 


再 记录 下 A10 手 环 采 集 得 到 


tH: 算法 在 慢 走 


两 者 进行 对 比 。 对 


时 平均 计 


手臂 | 


Ls 


度 来 看 ， 


为 了 解决 现 
不 同 用 户 行走 特征 的 问题 
的 计 步 算法 ， 在 数据 采集 阶段 采集 
率 行 走 的 数据 ， 


上 臂 慢 走 ( 将 


\ 态 


所 以 这 种 状 7 


快走 时 平均 计 步 精 
手 环 放 入 口袋 中 ) 


| 步 精度 比 小 米 2 精度 稍 低 一 点 。 一 般 手 环 都 会 戴 在 
精度 降低 是 可 以 接受 的 。 从 总 体 角 


本 算法 的 平均 计 步 精度 会 比 小 米 手 环 2 的 平均 计 步 


精度 要 高 。 


4 ”结束 语 


算法 在 经 


算法 对 不 同 运 动 状态 适 


四， 本 文 设计 了 一 


种 基 


应 性 差 和 无 法 适应 
于 自 适 应 阔 值 


了 不 同 


六 
外 


] 户 月 


日 三 种 步行 频 


实验 结果 表明 在 慢 走 、 快 走 或 跑步 三 种 不 


Hl 


量 的 数据 测试 和 不 断 地 改进 后 ， 
ij 频率 下 的 平均 计 


基于 自 适应 双 阔 值 的 计 步 算法 


上 精度 可 以 达到 91.88% 以 上 , 并且 
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统计 了 算法 对 不 同 用 户 的 


一 AN 


| 之 


HH 


纲 


了 走 特 征 
下 、 打 字 、 敲 桌子 等 ) 问 题 的 干扰 ， 完 成 了 预 
在 算法 中 仍然 手 
采集 得 到 的 数据 ， 
方向 是 继续 完善 算法 ， 
步 方 面 作出 更 精确 的 判定 ， 提 高 计 步 的 精度 。 


有 很 好 的 适应 性 ， 


有 效 地 排除 了 非 行走 因素 (起 立 、 
其 目标 。 不 过 
了 需要 加 强 之 处 ， 比 如 对 于 手 环 放 在 口袋 中 
算法 精度 还 需 继续 提高 等 。 下 一 步 的 研究 
运用 机 器 学 习 上 的 分 类 方法 在 检测 候 


表 4 提出 的 算法 与 小 米 手 环 的 对 比 结果 
Table 4 Comparison between proposed algorithm and xiaomi 
bracelet 
Ee BA = 

手 环 ”测试 项 Se 误差 步 数 ( 步 ) 准确 率 
A10 甩 臂 慢 走 200 165 -35 82.50% 
小 米 甩 臂 慢 走 200 177 -23 88.50% 
A10 甩 臂 慢 走 200 147 -53 73.50% 
小 米 甩 辟 慢 走 207 77 -130 37.20% 
A10 甩 辟 慢 走 200 153 -47 76.50% 
小 米 甩 臂 慢 走 200 104 -96 52.00% 
A10 甩 臂 慢 走 200 173 -27 86.50% 
小 米 甩 臂 慢 走 200 142 -58 71.00% 
Al10 电 臂 快走 200 186 -14 93.00% 
小 米 甩 臂 快走 200 160 -40 80.00% 
A10 甩 臂 快走 200 165 -35 82.50% 
小 米 甩 臂 快走 200 190 -10 95.00% 
A10 甩 臂 快走 200 170 -30 85.00% 
小 米 电 臂 快走 200 121 -79 60.50% 
Al10 电 臂 快走 200 132 -68 66.00% 
小 米 电 臂 快走 205 173 -32 84.39% 
A10 ”不 甩 辟 慢 走 。 200 186 -14 93.00% 
小 米 ”不 甩 辟 慢 走 ”200 201 1 99.50% 
A10 ”不 甩 辟 慢 走 ”200 168 -32 84.00% 
小 米 ”不 甩 辟 慢 走 ”207 213 6 97.10% 
A10 “不 甩 臂 慢 走 200 173 -27 86.50% 
小 米 ”不 甩 辟 慢 走 ”200 242 42 79.00% 
Al0 ”不 甩 辟 慢 走 ”200 168 -32 84.00% 
小 米 ”不 甩 臂 慢 走 ” 200 177 -23 88.50% 
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